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I. INTRODUCTION AND SUMMARY 



Graph theory and its applications play an important role in many areas of scientific research, from pure mathematics, 
to physics, statistics, biology and social sciences. In particular random graphs, that is graphs extracted with some 
probability from a suitable statistical ensemble, are interesting as a mean to implement the intrinsic complexity and/or 
chaotic nature of many physical, biological and social systems |l|. 

Among graphs, trees (that is graphs without loops) play a distinguished role: they retain a deep interest and wide 
applicability while still being amenable to detailed analytic studies. Random trees appears in several distinct contexts, 
like polymer physics, critical percolation j2i] a-nd two-dimensional quantum gravity (branched polymers Q). 

Generically, in a random graph the local coordination is itself a random variable taking values according to some 
probability distribution; while in the classical Erdos-Renyi theory of random graphs P| this distribution is Poissonian, 
in several examples of "experimental" complex networks it turns out to be "scale-free" , that is with a long power-law 
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tail [5j- This implies that on an infinite scale- free graph nodes with diverging coordination are rather frequent, causing 
some subtleties in the application of the law of large numbers. 

In this paper we concentrate our attention to homogeneous scale-free random trees subject only to the physically 
natural constraint of being embeddable in a finite dimensional Euclidean space. This class of random trees has been 
studied already in 0,0, although from a viewpoint different from ours: while in these works only some statistical 
averages relevant to characterize the geometry of these trees were analyzed, we study here in depth the probability 
distributions of the basic observables which specify the intrinsic geometry. Our results apply therefore to any single 
"generic" (infinite) scale-free random tree and not just to their statistical ensemble. 

The outline of this paper is as follows. First of all we recall some definitions of graph theory and we introduce the 
algorithm used to build homogeneous random trees, summarizing known results. The local connectivity dimension 
is then obtained from the scaling properties of the probability distributions for the volume and the surface of the 
spheres centered at the origin of the tree-constructing algorithm. Details of the scaled probabilities are studied and 
their asymptotic behavior is determined. Then the probability distributions for the graph-averaged volume and the 
surface are considered; the average connectivity dimension is extracted from their scaling properties. 

Our results for the growth properties of a random tree can be summarized in very concise heuristic formulas as 
follows: consider a sphere (intrinsically defined in terms of the chemical distance alone) of radius r centered on some 
node; for large r the volume v of the sphere and the coordination fluctuation Az within the sphere are random 
variables simply related as 

V ~ (Az)2r2 

This relation apply to all random trees with bounded growth rate (see below) regardless of whether Az has a finite 
limit as w ^ cx) or not. In the former case, as it happens for instance on trees with bounded coordination, one 
reads immediately out the connectivity dimension dc = 2. In the latter case, one need to estimate the asymptotic 
dependence of Az on v. For scale- free trees with tail exponent 2 < /3 < 3 (see below for the proper definitions), one 
finds Az ^ i;(3-/5)/(/9-i)^ yielding the /3— dependent connectivity dimension 



Suppose instead that one is interested in the graph-averaged volume v obtained by averaging over all locations of the 
center of the sphere (for a finite tree with large volume V). The random variable v is related to the coordination 
fluctuation Az over the whole tree just as before, that is 

V ~ (Azfr^ 

Now the behavior in r is always quadratic, regardless of the way the coordinations are distributed. Of course, in 
the scale-free trees with 2 < /3 < 3 one has Az ^ y(3-/3)/(/3-i)^ making it impossible to consider the standard 
thermodynamic limit V ^ oo. Nonetheless, we flnd that a well deflned limit exists for the "renormalized" volume 
yy{l3-^)/{P-'L) _ Therefore, in all cases we conclude for the average connectivity dimension 

4 = 2 

Finally, in the last section, the analysis of the probability distribution for the effective squared mass, defined through 
the Gaussian model on the tree, allows us to rigorously determine also the local spectral dimension dg. It fulfills 
the long-conjectured Q relation ds = 2dc/{dc + 1) already verified in through a general but completely different 
argument. 

II. RANDOM TREES 

A. Generalities about graphs 

A graph G is defined by a set of nodes (finite or countable), pairwise connected by a set of unordered links. If the set 
of nodes is finite, its cardinality will be denoted by \G\. The coordination number (or degree, or simply coordination) 
of a node x is the number of its nearest neighbors and it is denoted by z^. The intrinsic properties of the graph are 
determined only by the connections between the nodes; the metric is given by the so-called chemical distance: the 
distance c?(x, y) is equal to the number of links of the shortest connected path between x and y. This distance is used 
to deflne spheres; the "volume" Vr{x) is deflned as the number of nodes of the radius r sphere centered at node x, 



and the "surface" Sr (x) is the number of nodes of the corresponding spherical shell. Here we shall consider connected 
graphs in the limit of infinitely many nodes. 

Given a function F defined on the nodes, its average value F is defined as the infinite radius limit of the average 
over spheres; if such limit exists, it does not depend on the center of the spheres provided F is bounded from below 
and the graph has a bounded growth rate [oj, that is the surface vanishes with respect to the volume in the infinite 
size limit. 

The properties of graphs can be described by various parameters; here we consider the connectivity dimension jioj l 
and the spectral dimension The connectivity dimension describes how the volume of spheres scales with its 

radius r for r ^ oo; we can define a local connectivity dimension dc using the spheres centered on any given node x 

Vr{x) - r'^" (f) 

ft is not difficult to show that dc does not depend on x when the graph has a bounded growth rate. We can define 
also an average dimension dc, if the average volume is used (provided the limit defining it exists finite) 

- (2) 

These two connectivity dimensions usually coincide, but on strongly inhomogeneous graphs they can be different |12| |. 

The spectral dimension is related to long time properties of random walks on the graph, as well as to many other 
physical properties such as the infrared behavior of the Gaussian model • On a generic connected graph G this 
model is defined by assigning a real-valued random variable (j)^ to each node x G G with the Hamiltonian 



n = 



l^ioJ2'^l + lT.^'^--'^'y^^ (3) 



{x,y) 



where {x, y) denotes nearest neighbors pairs of nodes and /zq > is a free parameter (the squared mass in field- 
theoretic sense). The spectral dimension is determined by the infrared /kq ^ singularity of the diagonal element of 
the covariance 

Sir^gUD) ^ (4) 

where ((•)) denotes standard Gibbs expectation values weighted with exp(— 7i). One can show that dg does not depend 
on the choice of the node x [l3| . As for the connectivity dimension, it is possible to split the spectral dimension into 
a local one and an average one. The definition in eq. Q evidently corresponds to the local one; the average spectral 
dimension ds characterizes in the same way the singularity of the graph avera ge o f {{(px))- Again, the two dimensions 
usually coincide, but on strongly inhomogeneous graphs they can be different |l2l |. 

Suppose now that a certain statistical ensemble Q of (infinite) graphs is given. As a consequence, connectivity 
and spectral dimensions become in principle random variables with probability distributions derived from the "mi- 
croscopic" one specifying Q. Consider the connectivity dimension: if the graphs of G are rooted, we can use the root, 
denoted by o, as the center of the spheres those growth define dc; then the corresponding volumes Vr{o) are random 
variables distributed with some probability Pr[wj.(o) — v]- If this probability has a well defined scaling behavior for 
r —f oo, that is 

Pr[vr{o)=v] ^p(^) 

then the random variable Vr{o)r^'^ has a well defined probability distribution in the limit r — > oo and we may write 
again 



as in eq. £3), identifying d with a non-fluctuating local connectivity dimension dc- In other words dc is a property 
of a "generic" specimen of the ensemble, that is a a property of G itself. A similar argument applies to the average 
connectivity dimension dc and to the spectral dimensions using the random variable {{4>x)) ^'Hd its graph average. 

An alternative approach is used in ref. ,6j , where averages in the canonical or grand-canonical ensembles of certain 
random trees {i.e. graphs without loops) are studied rather than a single "generic" sample tree. Two parameters 
are introduced to describe the intrinsic geometry of these trees, also called "branched polymers": the Hausdorff 
dimensions dn and d^. The former is related to how the average two-point distance on graphs with V nodes scales 
with the size V: 



{d{x,y))v-V'/'' 



(5) 



The latter (called local Hausdorff dimension in [g ) is related to the behavior of the two-points correlation function 
gy (r) (proportional to the number of couples of nodes which are at distance r) for distances which are large, but 
much smaller than v^^'^": 

g'^\r)''^'^r''--\ 1 « r « F^/^" (6) 

It is claimed in that dn and differ on a certain class of "exotic" random trees characterized by unbounded local 
coordinations. 

Actually, it is not a priori obvious if and how parameters such as dn and d^ are connected to the connectivity 
dimensions dc and dc, although it is common lore to identify the Hausdorff dimension with the connectivity dimension 
when no distinction is made about local or average dimensions. One of the aims of the present work is just to fully 
answer this question in the case of homogeneous random trees, as we shall see below. 



B. Generalities about random trees 



Homogeneous random trees are built by the random independent extraction of every node's degree from a given 
distribution fz- This process can be formulated as a Galton- Watson branching process in particular, since we 
are interested in trees with bounded growth rate, we must consider the critical Galton- Watson case. The given 
coordination distribution must be properly normalized 

z 

and the average coordination must be equal to 2 because of the condition of bounded growth rate (see ref. 0) 

(z)=Ez/.-2 (8) 

z 

Clearly these sums have a finite number of addends if the coordination is bounded, while they become series if it is 
unbounded as in the scale-free case. The series, however, must be properly convergent in order equations and © 
to hold, so that must vanish faster than for z oo. Let us also introduce the probability generating function 
5(A) 

5(A) = = + /2A + + • • . (9) 



which enjoys the properties 



5(1) = 1 , 5'(1) = 1 (10) 



easily derived from those of ]z- 

A branching process effectively ends when a shell (that is a given generation of sibling branches) is made only by 
nodes with coordination 1, so that the next shell is empty. The trees produced in this way are all finite because the 
surviving probability after r generations vanishes for r ^ cxd in the critical Galton- Watson process Since we 

are interested in the thermodynamic limit we have two possibilities 0. The first one is an explicit preconditioning 
on non-extinction: this means that modified branching probabilities are used to avoid finite trees, while keeping 
unaltered the properties of the infinite trees [3| ■ This is achieved if the root coordination is chosen with probability 
jz while on every successive shell the coordination of the first node is extracted with a modified probability distribution 
fz = {z— l^)fz- The other possibility consists in calculating probabilities conditioned on the number of nodes V of 
the resulting trees, and then take the limit V 00. In this latter case the root coordination must be extracted with 
the refined probability 

K = nA, ^ = (11) 



since the root has as many branches as its coordination while every other node has one branch less than its coordination 
[see section 3.1 of ref. for details]. 



C. Trees with bounded coordination 



In [9j and | 16| the geometrical and spectral properties of bounded random trees were determined using non-extinction 
preconditioning. The surface and volume probability, in the large radius limit, are shown to be function only of the 
scaled variables s/r and v/r^ , respectively. Moreover, after the discussion of the auto-averaging property, also average 
surface and volume are fully analyzed. These results allow the rigorous determination of the connectivity dimension, 
both local and average: dc = dc ^ 2. The study of the Gaussian model on such trees leads also to the determination 
of the spectral dimension d^ = 4/3. Using a simple scaling hypothesis it is shown that the relation ds = 2dc/{dc + 1) 
should hold on a wider class of random trees. 

An important universality property is also found: all the average values and probability distributions in the large 
radius limit depend on the fz distribution only through its second moment or, equivalently, This parameter 

quantifies the coordination fluctuations, since 

z 

and characterizes all asymptotic behaviors. Altogether, these results can be summarized as 

v^g"il)r^ (12) 

meaning that v/[g" {l)r'^], as r — > oo, is a random variable with a well defined universal probability distribution. 

The approach based of the grand-canonical ensemble of branched polymers yields the result dn — dL — 2 when the 
coordination of each node is bounded. This corresponds to the "generic phase" of random trees according to @ . 



D. Scale-free trees 



In this paper we turn to scale-free trees, that is trees whose coordination distribution has a long power-law tail: 

fz — A z^^ , for z oo 

where (3 > 2. When the coordination is unbounded, the function g{\) may have singularities; in the scale-free case it 
becomes singular for A = 1. For any non-integer exponent^ its expansion contains a term proportional to (1 — A)''"^, 
that is: 

g{X)=g,{\)+C,{l-\f-^ (13) 

where 5a (A) is analytic in A = 1 

ga{\) = 1 - (1 - A) + C2(l + + Cfc(l - A)'= + . . . (14) 

and c, — AT{1 -~ (3). Therefore, if 2 < /3 < 3, g"{X) diverges for A ^ 1; since we found g"{l) to be an important 
parameter in the bounded coordination case, we expect this divergence to have many important consequences. As a 
first example it causes the divergence of {vr{o)), the expected volume of the balls around the root of the branching 
process. This can be intuitively explained, using non-extinction preconditioning, because on every shell the coordi- 
nation of one node is chosen with probability distribution fz, whose first moment is infinite. If auto-averaging holds, 
then also would diverge on a single scale-free tree and the definition ^ of average connectivity dimension appears 
troublesome for scale-free trees. Thus, in order to deal with finite graph-averaged quantities, we will also use the finite 
volume approach, regarding the volume ^ as a regulator. 

We may then find an heuristic argument based on eq. p2|l to determine the connectivity dimensions dc and dc for 
scale-free random trees. First of all, if /3 > 3 nothing changes because g"{l) is still finite, so that dc = dc = 2. On the 
contrary, if 2 < /? < 3, we have to face the divergence of 5"(1). But coordination fluctuations are finite for a finite 
number N of nodes, and diverge with iV in a way fixed by the long tail of fz, since all coordination extractions are 
independent. Let Zmax(-^) be the largest coordination extracted, estimated from 

oo 

V f - z^-^ = — 



^ for integer /3 the singularity has a different, logarithmic, form; however, since the final results will be analytic in /3, one is able to extend 
them also to integer values of /3. 



then we may estimate fluctuations by 



where we defined the exponent v as 

r.= (3-/3)/(/3-l) 
Thus for the local growth we have the consistency relation 



which implies 



V 1 



y ^ ^2/(1-1') = ^(/3-l)/(/3-2) 



P-2 



When the graph-averaged Vr is considered for trees of V nodes, the coordination fluctuation over the whole graph 
should be used instead: 

^^^2^(3-/5)/(/3-l) 

From this we read the scaling exponent 2 independently on V , but without the standard thermodynamic limit. Of 
course this is not a real proof and a more rigorous approach confirming these results is adopted in the following 
sections (see also Appendix 0; at any rate, these exponents are the same as those calculated for scale- free branched 
polymers with 2 < /3 < 3 in the grand-canonical approach: d/f = (/? — — 2) and = 2 |^. 

III. GROWTH STATISTICS 

A. Scaling probabilities for local surface and volume 

The probability Pr [sr(o) — s\ (Pr[i;,.(o) ~ v\) for the surface (volume) of the radius r sphere around the root of an 
infinite tree can be found, recursively on r, using the modified probabilities given by explicit preconditioning (see 
where a slightly different notation is used). The corresponding generating functions 

G^(A)=^A^Pr[s.(o)=.] 

S 

G;!(A) =^A" PrK(o) =w] 

V 

satisfy functional recursion rules easier to analyze, given by equations (15) and (16) in that is 

G;!+i(A) = A.g,+i(A).9;(A) 
G^+i(A, 1) fer+l(A) 

^^^-A5(/^.-i(A))^^ 

where 

gr+i{\) ^ g{9r{X)) go(A) = A (15) 

hr+i{X) = Xg{hA\)) ho{X) = X (16) 

These recursions involve only g{X) and g'{X), which are finite for < A < 1 for every value of /3 > 2, and therefore can 
be solved numerically with high accuracy for any r (subject only to computational limits), thus providing a solution 
of our problem from a practical point of view. In figure^the r — 10 volume probability corresponding to a tree with 
(3 = 2.5 is plotted against the "experimental" frequency distribution; this is obtained by sampling, over 1000 different 
graph of 8 • 10^ nodes each, the volume of 2000 balls centered on randomly chosen nodes. The agreement is good, 
and we notice the power law tail with exponent that can be estimated as ~ —1.5. In other cases (not shown) this is 
always equal to 1 — /3. 
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FIG. 1: Numerically calculated volume distribution vs. simulation histogram 

Analytically, upon substitution of eqs. (I13|l and (|14|) for g{X), one can solve the recursions rules H15|) . Hlt)|) for the 
first terms of the expansions of gr{X) and hr{X) as A ^ 1; 

griX) ~ 1 - (1 - A) + C2r(l - A)^ + . . . + c,r(l - Xf-^ 

hr{X) ^ 1 + r(l - A) + ^r3(l _ A)2 + . . . + ^r^(l - Xf-' ^^^^ 

where only the leading order term in r is shown for every power of (1 — A). Similarly for G'^(A) and G^{X) one finds 

G^(A) c 1 - 2c2r(l - A) + ... - (/3 - l)c,r(l - Xf^^ 
G^(A) ~ 1 - C2r2(l - A) + . . . - c,r'3-i(l - A)''-2 

Through the rules of (inverse) discrete Laplace transform, the singular terms with noninteger powers in 1 — A deter- 
minine power law tails in the probability distributions at fixed r: 

Pr[sr{o) = s] Ars-^0-''> , PrK(o) = v] ™ -^r^-\-(0-^) 

(3-1 

in agreement with numerical data. The local connectivity dimension dc is determined if we find scaling forms of these 
probabilities such that, as s — )■ oo or i; — > cxd and r ~+ cxd. 

This means that dc must be such that Gf.{exp{—u/r'^''~^) and G5!(exp(— ^/r'*'^) have finite non-trivial limits for r ^ oo. 
The right scaling exponent is found by examining equation p8(l : if /3 > 3 the leading order term in (1 — A) is the 
linear one so that finite functions of u and ^ are obtained with 

4 = 2 for /3 > 3 



while if 2 < /3 < 3 the most important term is the singular one so that the correct scaling is given by 

dc = for 2 < /3 < 3 (19) 

In appendix^we present all explicit calculations about the surface and volume probability distributions, for the two 
cases /3 > 3 and /? < 3, which confirm this simple scaling analysis. 

B. Fixed volume expectation values 

In this section we study the local growth properties of scale-free random trees of finite volume V . We shall 
reconstruct the dependence on V of the surface and volume probability distributions from the expectation values of 
Sr{o) and Vr{o) and their powers. We are lead to consider expectation values over finite-size trees, regarding their 
volume y as a regulator, because the power-law tails with exponent 1 — /? found in the previous section for Pr[sr(o) = s] 
and Pr [11^(0) = v\ imply the divergence with V of high enough moments of Sr(p) or Vr{o). 

Let us observe that in this context the average is made over all possible rooted trees with V nodes; this can also 
be seen as the graph average over all nodes of an unrooted tree, which is then further averaged over all possible 
realization of the unrooted tree, that is 

{0)v = {0)v (20) 

for every observable O. In particular 

{Sr{or)v ^{-^)V 

with a similar relation for the volumes. No dependence on the root o survives and we shall drop it from the expectation 
values. ^ 

Let us start by studying the probability P{V) that any rooted tree T with V nodes is produced by the branching 
process without the non-extinction precondition 

P{V) = Pr[|T| - V] 
By the tree-producing algorithm, P{y ) satisfies the identity 

^(^) = E/-E '^(^-i-E^.)flm) 

^ Vi,...,V, j=l j=l 

where P{V) = Pr[|i3| = V] is the same probability relative to a branch B, that is a tree whose root has an incoming 
link (and z — I branches) and whose coordination is extracted with probability ; it satisfies a similar equation with 
z — I instead of z and fz instead of fz : 

z Vu...,V,-i j = l j=l 

The generating functions G(A) = J2v P{V)X^ and G(A) = P(y)A^ satisfy 

G(A) =A5(G(A)) , G(A) = A5(G(A)) (21) 
where g{X) is the usual probability generating function defined in equation J^, while 

z 

Normalization of fz implies g{l) = 1 and from the definition easily follows that g{X) derivatives are proportional to 
those of g{X); more precisely we have 



5W(A) = 7Vig('=-i)(A) 



Since we are interested in the V oo limit, which corresponds to A ^ 1, we can use expansions (I13|) and (|14|) to 
obtain the asymptotic behavior of G{X) and G(A). For different values of /3 we have different situations depending on 
which term, the singular one or the quadratic one, is lower order. 
For /3 > 3 we have 

G(A)^f-c-^/^(l-A)V2 
G(A) ~ f- Aric^'/'(1- A)i/2 

so that, in the V oo, the probability P{V) reads 

P(V) ~ N, I y^/^ 

2^C27r 

For 2 < (3 < 3, instead, the singular term produces 

G(A)^l-c:^/(''-^)(l-A)i/('^-i) 
G(A)^f-7Vicri/(^-i)(l-A)i/(''-i) 

Thus we have, for V ^ oo 

Let us now turn to the calculation of the expected size of a shell of radius r in trees with volume V. First of all let 
us introduce the joint probability Pr{s, V) that the size of the tree is V and the r-th shell has s nodes: 

P,{s,V)=FT[sr{o)=s,\T\ = V] 
It is related to the same probability relative to branches Pr{s, V) as 

z z z 

Pr{s,V)^J2f^ E E S{V-l-J2V,)S{s~J2'^)I[Pr-l{s,,V,) 
z Vi,...,V^ si,...,s^ j = l j = l j=l 

while Pr{s, V) satisfies a similar recursion which is obtained from the previous one by replacing z with z — 1 and fz 
with fz, 

z-\ z-l z-1 

= J2 E S{V-l-Y,V,)Sis-Y,s,)l[Pr-i{s„V,) (23) 

z Vi,...,v;_i si,...,s2„i j=i j=i j=i 

Clearly the expectation value of Sr{o) on trees with V nodes is now given by 

is.)y = = SM (24) 

PiV) PiV) 

where Ey[s] is defined as the weighted sum in the first numerator. Analogously we can define i?Jf'[s] as 

S 

and the two generating functions 

V V 

A recursion equation for the latter can be found using equation H23|l 

i?,^M = ^A^^.p,,(.,y) 

V s 

= ^E/^ E E (E^^)n^"^^^-(^-^.) 

z yi,...,v'^_i fe=i j=i 

^\Y,Mz-i)E^^MGi\r-' 

z 

^Xg\G{X))E^,M 



and one for [s] is similarly obtained 



Since so{o) = 1, i?o'[s] is equal to P{V), so that the initial condition reads 



Then, for r > 0, the solutions read 



E^[s] = G{\) 



E^[s]^GiX) Xg'iG{X)) 



E^[s]^mXG{X)g{G{X)) A.g'(G(A)) 



(25) 



Now {sr)v can be easily found, but we have to consider the cases f3 > 3 and /3 < 3 separately. 
For /3 > 3 the A ^ 1 asymptotic expansion of E^ [s] reads 

E^[s] ~ iVi [l - (2 + 2c2(r - l))c2"'/'(l - X)'/^' 

which imply the large V limit of {sr)v 

{sr)v ~ 2 + 2c2(r- 1) 

For 2 < < 3 we have, instead, 
and 

{sr)v ^h{(3){r-l)V' 

where 



(26) 



k, (0) -(0- r(-i/(/3-i)) 

^^^^^"^^ r(-(/3-2)/(/3-l)) 

In the first case the V ^ oo limit is finite, while the second one we find a divergence, as expected. To understand 
how equation (|26() relates with the known dc, we now calculate higher moments of Sr- First of all we generalize the 
notation P,.(0,y), E}f[0] and E^[0] to every observable O 

Pr{0,V)=PY[Or{x)^0,\T\^V] 

Er[0]=J20PriO,V) 
o 

E^[0]=J2X''E^[0] 

V 

together with their branch counterparts Pr{0, V), Ejf[0] and E^[0]. Now it is easy to show that the recurrence rule 
for -E^ [5^] reads: 

E^[s'] = Xg\G{X))E^_,[s'] + Xg"{GiX)) {E^^.^f 
Thus for 2 < < 3 the second moment is found to read 

(s2),,~fc2(/3)(r-l)V"^+i/(^-i' 

where k2 (P) is some function of only. 

For all higher moments £"^[5"] we find, as a general rule, that the single E^_j^[s"] on the right hand side is multiplied 
by Ag'(G(A)) while the products of k are multiplied by Xg'^'^^ {G{X)). Moreover the equations for £';^[0] are the 



same as those for E^[0] with g{-) substituted by g{-). Therefore we obtain for aU higher moments, up to the leading 
terms in the A — s- 1 Umit, 

E^[s-] Ag'(G(A))i?^Js"] + Ag(")(G(A))i?^iM" + . . . 

This entails the final result 

~fc„(/3)(r-l)y''+("'i)/('5-i) 

which is compatible with a scaling form for the s,- probability just given by a finite-size scaling of the one found in 
section IlII Al that is 

where dc is the already known local connectivity dimension 

dc — — 

^ 13-2 

and the hmction q{x^ y) is such that 

'a;"g(a;,j/)da;~y('5-2-«)/(/3-i) 

Similarly it is possible to calculate the finite volume moments of the volume Vr{o) which read: 

«)y ~ r"+i T/("+2-«/(/3-i) (28) 
The proper scaling compatible with this result is 

Pr{-\y)^^j[-^^^^) (29) 

if the moments of the scaling function y) behave as 

a:>(x,7;)dx~2/('5"2-n)/(/3-i) 



Scaling H29|l agrees with numerical simulations as shown in figure |2] For several different values of r and V, we 
built 720 trees and calculated the volume v of radius r balls around 1000 (for each tree) randomly chosen nodes. Thus 
different "experimental" probability distributions are found; we considered the integrated probability Qr{v\V) = 
T^v'Kv Pr{v\V) and plotted the values of v corresponding to Qr{v\V) — 0.1, 0.2, 0.4, 0.6, 0.8, 0.9; (after the proper 
rescaling by r''-") versus r'^" /V . For everyone of the six values of Qr{v\V), the points corresponding to different values 
of r and V fall roughly on a curve, apart from statistical errors (fluctuations are expecially large for Qr{y\V) = 0.9, 
due to the long power-law tail of Pr{v\V)), so that the form l|29|l of the probability scaling is supported. 



C. Probability of the average surface 

In the bounded coordination case we were able to prove the auto-averaging property for every local observable: the 
graph average O is non-fluctuating on infinite trees and coincides with the expectation value (O) in the branching 
process around the origin 9, 16] . This implies immediately that local and average connectivity dimensions coincide. 

In the scale-free case the situation is more involved because many relevant expectation values do not even exist on 
infinite trees, so that we cannot prove the auto-averaging property in the same way. Thus, to determine the average 
connectivity dimension d^ we will try to directly investigate the probability distribution of the average surface s7 in 
the large r limit. In particular, we look for a scaling form of the type 
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FIG. 2: Group of symbols of the same color (online version) or shade of gray (printed version) correspond to all possible 
combinations of r = 10,20,30,40,50 (using the same symbol) and V = 1,2,4, 8 ■ lO'^ (using different symbols as explained in 
the legend) . 



A very close approach is possible also for the average volume Iv, but it is more involved and will not be reported here. 

The results of the previous section, together with equation (|20|l , allow us to calculate only the first moment of such 
probability distribution. We now reconstruct it (at least its scaling form) through the calculation of higher moments. 

This can be done by introducing Sr, defined as the sum of s,-(a;) over all nodes x of a tree: 

Then on a finite-sized tree, the average surface s^, is given by 





and its (finite volume) moments are simply obtained 



yn 



The recursive nature of the trees allows us again to find the composition rule for Sr in terms of the values of Srj 
and the local surface Skj (with k < r) relative to j— th branch around the root o. ft reads 



SrJ + 2Sr-l,j + ^ ^ Sk,jSr-k-2,j' 



(30) 



The recursion rule for the joint probabihty of Sr, si, . . . , Sr~i and V easily follows from (|30|l and the recursion rules 
for Sr and V. With calculations similar to those of the previous section one obtains 



r-2 



E^[S] = Xg'{G{X))E^[S] + 2Ag'(G(A))£;^i [s] + Ag"(G(A)) ^ i?,^[s]£;^,_2M 

k=Q 

and, as usual, the equation for E^[S] which is equal to this one with all g(-) substituted by g{-). The expectation 
value 

\Sr)v 



VP{V) 

is then obtained by straightforward calculations in the 1/ ^ oo limit, and the result is equal to the large V limit of 
{sr)v, as expected [recall eq. (jSOJ]- 

The differences of course emerge when higher moments are considered, since 7^ i'^r)" for n > 1. ((s7)^)y can be 
obtained using the equations (only the leading terms are displayed) 

E^[S^] c^Xg'{G{Xj)E^[S^] + ... 

E^[S^] ^Xg'{G{Xj)E^[S^] + X g" {G{X))E^[S]^+ 

■r-2 r-2 

+ Xg^^HG{X)) J2 E E^[s]E^,[s]E^_,_,[s]E^_,,_,[s] + ... 

k=0 k'=0 

and turns out to read 



mr)v 



(2 + 2c2(r- for/3>3 
k2iP)ir - l)^V^'' for2</3<3 



(31) 



where k^i^P) does not depend on r or V. 

Then all the following moments can be calculated through the equations 

E^[S^^]^Xg'{G{X))E^[S'^] + ... 

E^[S^^] ^ Xg'{G{X))E^[S^^] + Ag"(G(A))i ^ u ■., E^[S']E^[Sr-'] + 

+ A5"'(G(A))i E ,,,v,(,!'^._^v), ^.'[^^]^.'[^^'>.'[^"'^"'^1+ 
+ ... + Ag(")(G(A))S,^[^]"+ 

+ A5(^"nG(A)) E ^fe.W^'-fe.-2W---<W^^r-fc„-2W + --- 

All ,. . . , /cji 

From the leading singularity in 1 — A we read their asymptotic behavior for large V: 

^ /(2 + 2c2(r-l))" for/3>3 
^^ "^^ " \fc„(/3)(r - 1)"\^'- for 2 < /? < 3 

For /3 > 3, since all moments are powers of the same finite quantity, the probability distribution for must be a 
delta function in the y — > oo limit 

Pr[s7 = s] ^ S{s - (2 + 2c2(r - 1)) = S{s - (s,)) 

Thus in the thermodynamic limit the average does not fluctuate and coincides with the expectation value for the 
surface around the root of the tree and the auto- averaging property is verified for the observable Sr- In particular we 
can extract the average connectivity dimension 



Similarly the auto-averaging property can be verified for the volume; we expect it to hold in general for every observable 
with a finite expectation value. 

When 2 < /3 < 3 instead, from eq. (|32|l we can read the following scaling form for the average surface probability 

^^(^1^) - ^-^{^) (33) 

This shows that for every V there exists a scaling form of the probability distribution of s for all scale-free trees of 
size V, but it does not have a non-trivial standard thermodynamic limit. Although s7 diverges on infinite trees and 
the standard definition of average connectivity dimension cannot be applied, we see that a IZ-independent probability 
exists for the "renormalized" average surface sV~'^ and that this scales according to the average connectivity dimension 
dc = 2. 

For finite trees a comparison can be made with the definition of (equation ©): since the average surface differs 
from the two point function only through a normalization factor 

gy\r) oc {sr)v oc r 

our calculation not only confirms the value = 2; it also shows that this is not a property of the first moment only, 
but it holds "in probability" , that is for every "generic" scale-free random tree. 



IV. SPECTRAL DIMENSION 

In order to find the (local) spectral dimension of scale-free trees, we may repeat the approach of [T^ substituting 
fixed radius averages with fixed volume averages. The quantity ((0^)) is given by a normalized Gaussian integral over 
all variables y G G; if we perform all integrations except the one over (jj^ we are left with a last integral which is 
also Gaussian and normalized, thanks to the self-reproducing property of Gaussian integrals. Therefore we can define 
the effective squared mass ^J.{x) from the width of this last integral: 



2tt 



1 



On a tree produced by a branching process the rules of Gaussian integration allow us to express the effective squared 
mass of the root as a function of those of the branches [T^ . 



y branches 
of : 



(34) 



Similarly the effective squared mass of a branch can be expressed as a function of those if its sub-branches, in a 
recursive way. Since the recursion rule for ^j.{x) is highly nontrivial we are not able to find directly its probability 
distribution or its moments; instead we expand ^i{x) in powers of /io and we consider the coefficients Vn 



.(x) = £(-l)"+iK(a;)Mo" 



Substituting this expansion in equation H34|) . recursion rules for the cofficients are obtained: 



y branches 
of X 



V2{x)= ^ [V2{y) + VM' 



y branches 
of : 



and, in general 



Vn{x) 



y branches 
of X 



where Tn is explicitly given in appendix ^ The important feature which allows the calculation of Vn moments is 
that the expression of V^(a;) is linear in the 
Clearly V\ is just the volume V so that 



(^1 



and 



Now we need the averages for n > 2. We use the same notation as in the previous sections (P{^, 'E''^V\^ 

P(-), . . . ) with the only difference that now observables do not depend on r. 

Let us start with the expectation value of Vi\ equations for £''''[V2] and £''''[V2] can be easily obtained in the same 
way as before 



= \~\g\G{\)^'\g\G{\))E\vl\ 



(35) 



This is valid for any distribution, but now the 2 < /3 < 3 and /3 > 3 cases must be examined separately. 

When 2 < /3 < 3 we can read the asymptotic behavior of (V2) from the leading 1 — A singularity in equation H35|l : 



p{y) 



y(3/3-4)/(/3- 



In a similar way it can be proved that (see appendix IbI for details): 



yl+(n-l)(2/3-3)/(/3-l) 



y -> CX) 



We can then write 



(36) 



(37) 



The existence of the thermodynamic limit [l^ requires that F\(t) ~ i-(/3-i)/(2/3-3) ^^j. ^ gQ \\^qx the powers of 

V cancel out and a finite limit is obtained. After the y — > 00 limit, for /ip ^ we have (see section 5 of 16] for a 
discussion on the order of the limits) 

(/^)oo = lim {ii)v oc ^ /io ^ 

Equation ljB3|l in appendix ^ can be used to show that the higher moments of [i follow similar scaling laws: 



so that 



n(/3-2)/(2/3-3) 

Mo 



Mo ^ 



Since all moments are proportional to powers of the same quantity, we can say that there exists a limit probability 

-3) 
' SO 



distribution for the scaled variable /i/Xp 2)/(2/3 3) \\^qx the local spectral dimension reads: 



'2/3-3 



for 2 < /3 < 3 



When /? > 3 equation allows us to write 



(■t^2)y = 



E''\V2\ 
P{V) 



y5/2 ^ 



V ^00 



All other moments are calculated in appendix iBl and, following the same steps as before, we can write 



so that 



ds = ^ for /3 > 3 



as in the bounded coordination case|l6j. 

Finally we notice that in both cases the relation 



„ dr 
4 - 2- 



dc + 1 

between the connectivity and the spectral dimension is fulfilled. 



APPENDIX A: DETAILS ON PROBABILITY DISTRIBUTIONS 



In order to find the probability distributions for surface and volume we can follow the same steps of 9]. If the 
proper scaling form of the surface and volume generating functions are substituted in their recursion rules, differential 
equations are obtained by consistency requirements. The solutions are the Laplace transform of the scaled surface 
and volume probability distributions. For /3 > 3 both Laplace transforms can be calculated, while for 2 < /3 < 3 only 
the surface probability Laplace transform and the asymptotic behavior of the volume one can be obtained. 



1. /3 > 3 



When /3 > 3, dc is equal to 2 and we can write gr{X) as 

g,(e-"/'-) =1-^ 

with the condition 



a,(u) 



r/3-2 



a'(0) = 1 



a(0) = , 

that follows from gr{l) — g'j.{l) ~ 1 and 
derived from equation l)17|l. This expansion imply those of G^(A) and Pr[sr(o) = s] 



for 



(Al) 

(A2) 
(A3) 



G^(e-"/'^) =a'{u) + ...- 
1 



Pr[sr(o) = s] = 



1 



(A4) 



where ^(a;) and are the inverse Laplace transforms of a'(u) and a'^(u), respectively. The analytic term (f>{x) 

is the same as in the bounded case; the singular term 4'*{^)i even if it is suppressed by powers of r, is essential 
in order to reproduce the power-law tail. The functions a{u) and a,(w) are determined by inserting equation (|A1|) 
in the recurrence H15|l . and equating the coefficients of the corresponding powers of r. Thus two coupled differential 
equations are obtained (see ref. for details on the method in the case of bounded trees, when only a{u) is required): 

ua'{u) ~ a{u) — ciaiuf' 

ua'^{u) = [(3 — 2)a*(u) — 2c2a{u)a^{u) + c^a{u)^^^ 



The solution of the former, obtained with the use of the initial condition (jA2|l . reads 



a{u) = 



1 + C2U 



while the general solution of the latter reads 



a*(u) 



C2(/3-4) 



1 



(1 + C2uf (1 + C2Uy^-^ 



The value of k is found by comparison with the condition ljA3|l and turns out to be /c = 1. The inverse Laplace 
transforms (j){x) of a'{u) is 



0(x) 



-X/C2 



while only the asymptotic a; — > oo behavior of (j)^{x) can be explicitly calculated: 

Therefore, even if the contribution of </>,(a:) in equation IIA4p is suppressed by powers of r, it is the most important 
one for large x because of its power- law tail compared with the exponential one of (j){x). Moreover, if /3 > 4 there 
are other lower order corrections before the singular one; however since their contribution in the Laplace transform 
G^(exp(u/r)) is analytic for m — > 0, they are exponentially vanishing for large s/r. The large s tail is therefore due to 
the singular term. 

Similar calculations can also be done for the volume probability. The first step consists in writing the expansion of 
hr{\) with the proper scaling: 



^(3-2 



with the condition 



6(0) = , 6'(0) = 1 



that comes from ^.^(1) — 1 and ^.'^(1) — r and 



(A5) 



(A6) 



(A7) 



derived from equation H17|l. Inserting equation (|A5|) in the recurrence for hr{X) (equation (|16|l ) and equating the 
coefficients of corresponding powers of r, we obtain: 

2^6'(0-6(e)-C26(C)'+e 

= (/? - 2)6,(0 - 2c26(06,(0 + c,6(0''-' 

Their solutions, imposing the conditions (|A6p and HA7|I . can be found to read: 

T 



&*(0 



tanh y' 



2i 



tanh ■\/c2^ 



\ _ /? - 2 2^1 (1, /3/2; 1 + /3/2; (tanh ^f) 
13 (coshV^)2 



where 2-^1(0, 6; c; z) is the hyper-geometric function; in our case it can also be written as: 



2Fi(l,6;l + 6;z) = ;^ 



(A8) 



71=0 



The volume probability Pr[i;j.(o) = v] and its generating function G'JJ(A) now can be written as 

G:{e-^n = m (1 + ...-^ + ...) 



Pr[vr{o) = u] = — $( — 



r/3-1 



$4^ 



where $(x) and $,(a;) are the inverse Laplace transforms of and ^(0'*(0; respectively. Probability normalization 
requires l{0) = 1 and U{Q) = 0; the functions /(O and /,(0 related to 6(0 and 6*(0 by the differential equations: 

^no = -c2fc(o^(«) 

2^^:(0 = (/? - 3)/*(0 - 2c26,(0 + c,(/3 - 1)6(0^"' 



The solution of the former is found in |9( and reads 

m = 

while for the latter we obtain 



(cosh y/c2S,Y 



(/3-3)/2 



P-l /3 + 1 



= - - 

C2 \C2 



where t = tanh y/c^. The solution (^) can also be written as a power series using equation (jA8|l 



(/3-3)/2 



^^+i+2»(2n + 2) 
/? - 1 + 2n ' ^ (/3 + 2n)(/3 + l + 27i) 



C2 VC2 

from which we can extract the asymptotic behavior for vanishing ^ 

These results imply $(a;) has an exponentially vanishing a; ^ oo tail (see 0) 



$(2;) ~ —5- e 



>^/(4c2) 



while $*(a;) has a power-law one 



Therefore we may repeat the same consideration as in the surface probability calculation: even if <i>*(x) is suppressed 
by powers of r, it is the most important term for large v/r^ . 



2. 2 < /3 < 3 



In this case the expansion of gr{^) reads 



/(/3-2) 



a{u) 



with the conditions 



i(0) = and a'(0) = 1 



(A9) 



(AlO) 



The surface probability distribution and its generating function G^(A) can then be written as functions of scaled 
variables as 

G^(exp(-u/»-''^-i)) = a'(u) + ... 

where is the inverse Laplace transform of a!{x). As before, by substituting equation 1A9|I in the recurrence 
equation for (?r(A) and by equating the numerators of the subleading terms, we obtain the differential equation 

uo!{u) = a{u) -c4/3- 2)a{uf-^ 



The solution, using the conditions IjAlOp . reads 

a{u) - 

Now, from the asymptotic behavior of a'{u) 
a'{u) = \l + c,{l3 - 2)u'^-^ 



[1 + c4/3-2)m'3-2]i/(/3-2) 



(All) 



ai-i)/(/3-2) I 1 _ - l)u'^~^ for M ^ 

" 1 [c,(/3 - 2)]-('3-i)/('3-2)u-'5+i for u^oo 



we can determine that of the scahng function 



0(a:) 



c for a; ^ oc 

v^~^ for a; ^ 



which has the expected power-law behavior. 

Let us now turn to the volume probability. The conditions hr{l) — 1 and h'(r) — r lead to the the expansion 

m 



with 



6(0) ^ , 5'(0) = 1 (A12) 
The differential equation for 6(^) is found by substituting this expansion in the recurrence Ijlfill : it reads 

if3 - m'iO = 6(0 - c.(/3 - 2)b{0P-' + (/? - 2)e (A13) 
By direct substitution one can easily verify that 

6(0 = (e/c.)^/(^^^) 

is a solution, but does not satisfy the second condition of equation ljA12(l . However we obtain a separable differential 
equation if we let 6(0 = ^/^''"^^KO 



The solution is now implicitly given by 



oo ^ 

^ l + fc(/3-l) * 



fe=0 



This allows us to determine the asymptotic behavior of the solution 6(0 



e - (c* 7/3)0'- 1 forC-O 
for ^ ^ CX3 



The volume probability and its scaled generating function can now be written as 

G:!(exp(-eA'^^)) = ^(0 + --- 
1 / V 



/ V \ 



where is the Laplace transform of $(a;). As before Z(0) = 1 by normalization and Z(^) is obtained from the 
differential equation 

^l'{0 = -cJ3~2)l{0b{0^-' (A14) 
The small and large ^ asymptotic behavior of Z(^) can be obtained from those of and read 

forego 

- \ exp(-c,(/3 - for ^ ^ 00 ^ ' 



Therefore the asymptotic behavior of ^{x) reads 

$(a;) 

with 



X for X — > CX3 

a;~^/^ exp(— /ca;^"'') for a; ^ 



The probability distribution is exponentially vanishing for small volumes and it has a power law tail with exponent 
1 — /?, as expected. 

APPENDIX B: K MOMENTS 

First of all let us write the explicit form of J^niVi, V2, . . . , V^-i) 

n n — 1 n—1 k k 

fe=2ni = l nk = l i=l i=l 

Now we can write an expression for i?^[y„] and its generating function iJ^[V^] involving only Vi,V2, ■ ■ ■ , Vn-i- 

= E E*^!^! - 1 - E ^1.^) E {^''"^ [^"] + ^"""^ [^n(V^l, V-^, . . . , V;_i)]} n ^(Vi,/) 
z 

= Xg'iG{\))[E^[Vn] + E^[J^r.{Vi,V2, 

= "1 - Aff'(G(A))] '\g'{GiX))E^[J'„{VuV2, . . . ,K-i)] 

This shows that is more singular than E^[Tn{Vi, V2, ■ ■ ■ , Vn-i)], which has now to be calculated. To do this 

we need an expression for the generic multiple moment E^[Vi^V2 '^ ■ ■ ■ V^"-]. First of all notice that powers of Vi can 
be extracted, that is 

E^[V^Wt ■■■V^"] = V'''E^[Vt ■ ■ ■ Vj^"] 
The recurrence rule for E^[V^' ■ ■ ■ V^^] reads 

E^Vt ■ ■ ■ = [1 - Aff'(G(A))] {a5'(G(A)) ^ E'[V^^ . . . k,V^^-':F, ■ ■ ■ Vt\+ 

A5"(G(A)) E\vt'---v::-'']E^[vt-'---v::'^'^]+ 

(Bl) 



A5W(G(A)) flE\vt^---V::-'^\] 



where N = kj, the sums over kj^i are such that kj^i = kj, and only most singular terms are kept. Equation ljBl|) 
can be recursively used; the moments on the right-hand side either involve products with fewer terms or lower j V^ 's 
so that we end with moments of Vi, which are trivial. 

The general expression of E^[V^^ ■ ■ ■ V^"] for both 2</3<3and/3>3 cases 

^ [^1 ^2 •••v; ]^r^_^^._.j:^,^.(3,_i) (B2) 



can now be proved by induction. Substituting equation ljB2l) in the right-hand side of equation (|B1|I we see that all 
terms have the same singular behavior when 2 < /? < 3, while the g'(G(A)) and g"{G{\)) terms are the leading order 
ones for /3 > 3. 

The branch probability moments just calculated are related to the tree probability ones by 

E\Vn] = Xg'{G{X)) [E^[Vn] + E\Tn{VuV2. . . . , K-l)] 

^Xg'{G{\))E^[Vn] + ... 

and 

E^iV^^'' • • • V^>] = Xg'{G{X))E^[V^^'' ■■■Vj;"] + ... 

Therefore E-^[Vi^ ■••Kf"] has the same asymptotic behavior (equation HB2|) as E^[Vi^ This imply the 

following asymptotic behavior of E^[V^W^^ ■ ■ ■ Vj^"] 



'yEj fcj((2/3-3b-/3+2)-/3]/(/3-l) 2 < ^ < 3 



and finally 



'YT.j k,({2l3-i)]-l3+2)/(l3-l) 2 < /? < 3 

YkT.,k,(3j-i) ^^3 



{V^^v^^^ ■ ■ ■ - <! Ia y .".,(3,-1) : : o " ^ (B3) 
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